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Abstract 

Deeply virtual Compton scattering, 7*p --fp, has been measured in e~^p col- 
lisions at HERA with the ZEUS detector using an integrated luminosity of 
61.1 pb~^. Cross sections are presented as a function of the photon virtual- 
ity, Q^, and photon-proton centre-of-mass energy, W, for a wide region of the 
phase space, > 1.5 GeV^ and 40 < < 170 GeV. A subsample of events 
in which the scattered proton is measured in the leading proton spectrometer, 
corresponding to an integrated luminosity of 31.3 pb~^, is used for the first direct 
measurement of the differential cross section as a function of t, where t is the 
square of the four-momentum transfer at the proton vertex. 
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1 Introduction 



This paper presents cross-section measurements for the exclusive production of a real 
photon in diffractive e'^'^ interactions, ep — > e'-fp, a process known as deeply virtual 
Compton scattering (DVCS). In perturbative QCD, this process is described by the ex- 
change of two partons, with different longitudinal and transverse momenta in a colourless 
configuration. At the 7*p centre-of-mass energies, W, available for ep collisions at the 
HERA collider, for large momentum-transfer squared at the lepton vertex, Q^, the DVCS 
process is dominated by two-gluon exchange. Measurements of the DVCS cross section 
provide constraints on the generalised parton distributions (GPDs) [1-5], which carry 
information about the wave function of the proton) [2]. The transverse distribution of 
partons in the proton, which is not accessible via the F2 proton structure function, is 
accounted for in the dependence of the GPDs on the four-momentum transfer squared 
at the proton vertex, t. The initial and final states of the DVCS process are identical to 
those of the purely electromagnetic Bethe-Heitler (BH) process. The interference between 
these two processes in principle provides information about the real and imaginary parts 
of the QCD scattering amplitude [6-8]. However, the interference is expected to be small 
in the kinematic region studied in this paper [6,7]. 

The simplicity of the final state and the absence of complications due to hadronisation 
mean that the QCD predictions for DVCS are expected to be more reliable than for 
exclusive vector meson production which has been extensively studied in ep collisions at 
HERA [9-16]. Several measurements of DVCS at high W are available [17-20]. The 
analysis presented here is based on data in the kinematic range of 1.5 < < 100 GeV^ 
and 40 < < 170 GeV, an extension compared to the previous ZEUS measurement [18]. 
A subsample of the data in which the scattered proton is measured in the ZEUS leading 
proton spectrometer (LPS) [21] is used for the direct measurement of the t dependence of 
the DVCS cross section. 



2 Experimental set-up 

The data used for this measurement were taken with the ZEUS detector at the HERA ep 
collider in the years 1999 and 2000, when HERA collided positrons of energy 27.5 GeV 
with protons of energy 920 GeV, and correspond to an integrated luminosity of 61.1 pb~^. 
The subsample used to measure the t distribution was collected in 2000 and corresponds 
to an integrated luminosity of 31.3 pb~^. 



Hereafter, the positron is referred to with the same symbol, e, as the electron. 
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A detailed description of the ZEUS detector can be found elsewhere [23,24]. A brief 
outline of the components most relevant for this analysis is given below. 

Charged particles were tracked in the CTD [25]. The CTD operated in a magnetic field 
of 1.43 T provided by a thin solenoid. It consisted of 72 cylindrical drift-chamber layers, 
organised in nine superlayers covering the polar-angle [§ region 15° < 6 < 164°. The 
transverse-momentum resolution for full-length tracks was a{pT)/pT = 0.0058]?r©0.0065© 
0.0014/pr, with pt in GeV. 

The uranium-scintillator calorimeter (CAL) [26] covered 99.7% of the total solid angle and 
consisted of three parts: the forward (FCAL), the barrel (BCAL) and the rear (RCAL) 
calorimeters. Each part was subdivided transversely into towers and longitudinally into 
one electromagnetic section (EMC) and either one (in RCAL) or two (in BCAL and 
FCAL) hadronic sections (HAC). The CAL energy resolutions, as measured under test- 
beam conditions, were a{E)/E = 0.18 /Ve for positrons and a{E)/E = 0.35 /Ve for 
hadrons, with E in GeV. 

The position of positrons scattered at small angles to the positron-beam direction was de- 
termined combining the information from the CAL, the small-angle rear tracking detector 
(SRTD) and the hadron-electron separator (HES) [27,28]. 

The FPC [29] was used to measure the energy of particles in the pseudorapidity range 
rj ^ 4.0 — 5.0. It consisted of a lead-scintillator sandwich calorimeter installed in the 
20 X 20 cm^ beam hole of the FCAL. The energy resolution for electrons as measured in a 
test beam, was a{E) /E = (0.41 ± 0.02) /VE © 0.062 ± 0.002, with E in GeV. The energy 
resolution for pions was a{E) /E = (0.65 ± 0.02) /^/E © 0.06 ± 0.01, with E in GeV, after 
having combined the information from FPC and FCAL. The e/h ratio was close to unity. 

The LPS [21] detected positively charged particles scattered at small angles and carrying 
a substantial fraction, xl, of the incoming proton momentum; these particles remained 
in the beam-pipe and their trajectory was measured by a system of silicon microstrip 
detectors that was inserted very close (typically within a distance of a few mm) to the 
proton beam. The detectors were grouped in six stations, SI to S6, placed along the beam- 
line in the direction of the proton beam, between 23.8 m and 90.0 m from the interaction 
point. The particle deflections induced by the magnets of the proton beam-line allowed a 
momentum analysis of the scattered proton. For the present measurements, only stations 
S4, S5 and S6 were used. The resolutions were about 0.5% on the longitudinal momentum 
fraction and about 5 MeV on the transverse momentum. The LPS acceptance [22] was 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the proton 
direction, referred to as the "forward direction" , and the X axis pointing left towards the centre of 
HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity is defined as 
r/ = — In (tan |), where the polar angle is measured with respect to the proton beam direction. 
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approximately 2% and XL-independent for xl>0.98; it increased smoothly to about 10% 
as xl decreased to 0.9. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e'yp, 
where the photon was measured in a lead-scintillator calorimeter [30] placed in the HERA 
tunnel ai Z = —107 m. 

3 Monte Carlo simulations 

The acceptance and the detector response were determined using Monte Carlo (MC) simu- 
lations. The detector was simulated in detail using a program based on Geant 3.13 [31]. 
All of the simulated events were processed through the same reconstruction and analysis 
chain as the data. 

A MC generator, GenDVCS [32], based on a model by Frankfurt, Freund and Strikman 
(FFS) [33], was used to simulate the elastic DVCS process as described in [18]. The 
ALLM97 [34] parameterisation of the F2 proton structure function of the proton was used 
as input. The t dependence was assumed to be exponential with a slope parameter b set 
to 4.5 GeV~^, independent of W and Q"^. 

The elastic, ep — >• ejp, and quasi-elastic ep e^yY BH processes, where Y is a low- 
mass state, and the exclusive dilepton production, ep —>■ ee^e~p, were simulated using 
the Grape-Compto]S§ [35] and the Grape-Dilepton [35] generators. These two MC 
programs are based on the automatic system Grace [36] for calculating Feynman dia- 
grams. A possible contribution from vector meson electroproduction was simulated with 
the Zeusvm generator [37]. To account for electroweak radiative effects, all the generators 
were interfaced to Heracles 4.6 [38]. 

4 Kinematic variables and event selection 

The process ep e'-fp is parametrised by the following variables: 

• Q"^ = —q^ = —{k — k'Y, the negative four-momentum squared of the virtual photon, 
where k {k') is the four-momentum of the incident (scattered) positron; 

• = {q+pY, the squared centre-of-mass energy of the photon-proton system, where 
p is the four-momentum of the incident proton; 

Hereafter, the Grape-Compton generator is referred to as Grape. 
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• X = Q^/ (2P ■ q), the fraction of the proton momentum carried by the quark struck by 
the virtual photon in the infinite- momentum frame (the Bjorken variable); 

• Xl = the fractional momentum of the outgoing proton, where p' is the four- 
momentum of the scattered proton; 

• t = {p — p')"^ , the squared four- momentum transfer at the proton vertex. 

For the range of this analysis, > 1.5 GeV^, and at small values of t, the signature 
of elastic DVCS and BH events consists of a scattered positron, a photon and a scattered 
proton. The scattered proton remains in the beam-pipe where, for a subsample of events, 
it is detected in the LPS (LPS sample). 

The events were selected online via a three-level trigger system [23, 39]. The trigger 
required events with two isolated electromagnetic (EM) clusters with energy greater than 
2 GeV. The trigger efficiency was studied as a function of the lowest energy cluster, it was 
found to increase from 80% to 100% for increasing cluster energy and the Monte Carlo 
was reweighed according. 

The offline selection followed the strategy described in [18]. Two EM clusters were found 
by a dedicated, neural-network based, positron finder [40]. They were ordered in polar- 
angle and are in the following denoted as EMI and EM2, with di > 62- The first cluster 
was required to be in the RCAL with energy Ei > 10 GeV; the second cluster had to 
have a polar angle 62 < 2.85 rad and was required to be either in the RCAL, with energy 
E2 > 3 GeV, or in the BCAL, with energy E2 > 2.5 GeV. The angular range of the second 
cluster corresponds to the region of high reconstruction efficiency for tracks in the CTD. 
The association of a track discriminates between positron and photon induced clusters. 
For events with one track, a match was required between the track and one of the two EM 
clusters. Events with more than one track were rejected. To ensure full containement of 
the electromagnetic shower, the impact position of each EM cluster on the face of RCAL 
was required to be outside a rectangular area of 26 x 16 cm^ around the beam-pipe. 

The condition 40 < E — Pz < 70 GeV was imposed, with E = Ei + E2 and Pz = 
£'icos6'i + i?2Cos^2- This requirement rejected photoproduction events and also events 
in which a hard photon was radiated from the incoming positron. 

Events with CAL energy deposits not associated with the two EM clusters were rejected 
if their energy was above the noise level in the CAL [41]. In addition, the total energies 
measured in the FPC and in the FCAL were each required to be below 1 GeV [10,41]. 
These elasticity requirements also suppressed DVCS events and inelastic BH events in 
which the proton dissociates into a high-mass hadronic system. The sample was still 
contaminated by events in which a forward, low-mass hadronic system was not visible 
in the main detector. Alternatively, a clean sample of elastic DVCS and BH events was 
obtained by additionally requiring the proton to be detected in the LPS. 
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The LPS event was rejected if, at any point, the distance of the proton track candidate to 
the beam-pipe was less than 0.04 cm. It was also rejected if the X position of the track 
impact point at station S4 was smaller than —3.3 cm. These cuts reduced the sensitivity 
of the acceptance on the uncertainty in the position of the beam-pipe apertures. To 
suppress background from overlays of ep collisions with protons originating from the beam- 
halo, it was required that {E + Pz) + < 1865 GeV, where is the longitudinal 
momentum of the scattered proton. The variable was required to be within the range 
0.96 < xl < 1.02 to exclude non-elastic events [42,43]. The variable t was required to 
be in the range 0.08 < |t| < 0.53 GeV^ where the LPS acceptance was high and slowly 
changing. 

The kinematic region was 40 < < 170 GeV and 1.5 < < 100 GeV^. For the purposes 
of this analysis, the values of and W were determined for each event, independently of 
its topology, under the assumption that the EMI cluster is the scattered positron. This 
assumption is always valid for DVCS events for the range considered here. The electron 
method [44] was used to determine and the double-angle method [44] to determine W. 

5 Background study and signal extraction 

The selected events were subdivided into three samples, 

• 7 sample: EM2, with no track pointing to it, is taken to be the photon and EMI is 
assumed to be the scattered positron. Both BH and DVCS processes contribute to this 
topology. The sample consisted of 7618 events and 55 events after the LPS selection. 

• e sample: EM2, with a positive-charge track pointing to it, is assumed to be the 
scattered positron and EMI is the photon. The sample is dominated by BH events. 
The number of DVCS events is predicted to be negligible due to the large implied 
by the large positron scattering angle. This sample consisted of 11988 events and 33 
events after the LPS selection. 

• negative-charge-e sample: EM2, with a negative-charge track pointing to it, may have 
originated from an e^e~ final state accompanying the scattered positron, where one 
of the positrons escaped detection. This sample is dominated by non-resonant e~^e~ 
production and by J/ip production with subsequent decay into e"*"e~ and was used to 
study these background sources. It consisted of 764 events and only one event after 
the LPS selection. The diffractive electroproduction of p, uj and (f) mesons was found 
to be negligible [18]. 

In the kinematic region of this analysis, the contribution of the interference term between 
the DVCS and BH amplitudes is very small when the cross section is integrated over the 
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angle between the positron and proton scattering planes [6,7]. Thus the cross section for 
exclusive production of real photons was treated as a simple sum over the contributions 
from the DVCS and BH processes. The DVCS cross section was determined by subtracting 
the latter. 

The size of the BH contribution to be subtracted was determined using the e sample 
which consists of elastic and inelastic BH events and a small fraction of exclusive e~^e~ 
production. The exclusive e~^e" contribution was estimated with the negative-charge-e 
sample to be (6.4 ± 0.2)% and subtracted from the e sample. 

The inelastic fraction of the BH events was estimated from the difference in the azimuthal 
angles, A0, between the two electromagnetic clusters in the e sample. It was determined 
to be (16 ± 1)% and was negligible in the LPS tagged subsample [45]. 

The measured cross section of the BH process was (4± 1)% smaller than the expectations 
of the Grape program (a detailed discussion can be found elsewhere [45,46]). The Grape 
cross section was modified accordingly. 

The BH contribution to the 7 sample was determined by Grape and found to be (56±1)% 
for the untagged and (21 ±3)% for the LPS tagged sample. The BH-subtracted 7 sample 
was further scaled by (1 — fp-diss), where /p_diss is the fraction of DVCS events in which 
the proton dissociated into a low-mass state. Its value was taken, as in [18], from previous 
pubhcation [47], /p_diss = 17.5 ± 1.3+;^;^%. 

The W and distributions in the untagged sample (inclusive sample) and the xl and t 
distributions in the LPS sample, separately for the e sample, for the 7 sample and for the 
7 sample after BH and proton dissociation background subtraction, are shown in Fig. [H 
Also shown in the figure are MC expectations which describe the data well. 

6 Systematic uncertainties 

The uncertainties due to the reconstruction of the scattered positron and to the back- 
ground subtraction were evaluated by varying the selection criteria as follows: 

• varying the electromagnetic energy scale by ±2%; 

• restricting the E - Pz cut to 45 < E - < 65 GeV; 

• shifting the reconstructed position of the positron with respect to the MC by ±lmm; 

• changing the elasticity requirements by ±30 MeV in the EMC and ±50 MeV in the 
HAC sections; 

• changing the photon candidate energy by ±10%; 
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• varying the inelastic BH fraction by ±1%. 

Each individual systematic uncertainty affects single bins in and W typically by less 
than 5% and by less than 10% in all cases bar the highest bin where statistical fluc- 
tuations dominate. 

To evaluate the uncertainties due to the reconstruction of the final-state proton, 

• the cut on the minimum distance to the beam-pipe was increased to 0.1 cm; 

• the t range was tightened to 0.1 < \t\ < 0.5 GeV^; 

• the kinematic limit of the beam-halo background cut was lowered to 1855 GeV; 

• the X position of the track impact point at station S4 was restricted to —32 mm. 

The total systematic uncertainty was obtained by adding in quadrature the individual 
contributions. It was found to be ±8% on average which is smaller than the statistical 
uncertainties. 

For the inclusive sample the uncertainty in the determination of the integrated luminosity 
of ±2.25% and on the proton-dissociative background of are not included in the 

figures and in the tables. 

For the LPS data, there is an overall uncertainty of ±7% which originates mostly from 
the uncertainty on the simulation of the proton-beam optics. It can be treated as a 
normalisation uncertainty as it is largely independent of the kinematic variables and is 
not included in the figures and in the tables. It also includes the uncertainty on the 
integrated luminosity for the LPS sample of ±2.25%. 



7 Cross section determination and results 

The 7*p cross section of the DVCS process was evaluated as a function of W, using 
the expression 

where N°^^ is the total number of data events in the 7 sample in bin i of W and Q^, 
denotes the number of elastic and inelastic BH events in the 7 sample in the bin, and Nf^^ 
is the number of events expected in the 7 sample from GenDVCS for the luminosity of 
the data. The cross section as predicted by the FFS model is denoted a^^^{'y*p — > •yp) 
and was evaluated at the centre {Wi, Qf) of each and W bin. The differential cross 
section as a function of t was calculated from the LPS tagged sample for which /p_diss is 
zero. All the results are listed in Tables [T]-|H 
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The 7*p DVCS cross section, ct'^'p-^tp, is presented in Fig. [2] as a function of at 
W = 104 GeV and as a function of W at = 3.2 GeV^. The cross section shows 
a fast decrease with Q^. A fit to the dependence of the cross section, assuming 
the functional form a'^ p^^p[Q'^) ~ Q^^", was performed for W = 104 GeV yielding 
n = 1.54 ± 0.05(stat.), smaller than expected for a pure propagator term [33]. The re- 
sult is in agreement with other DVCS measurements at HERA at lower W [17-20]. As 
expected for DVCS [33], the decrease of the cross section with is slower than for 
exclusive vector meson production [12-16,48]. The cross section increases with W. In 
pQCD-based models, this behaviour is related to the increase of the gluon content of the 
proton with decreasing Bjorken-x. A fit to the W dependence of the cross section, as- 
suming a functional form a'^*P^'^P(W) ~ W^, was performed for = 3.2 GeV^, yielding 
6 = 0.52±0.09(stat.). This result is in agreement with the previous measurements [17-20] 
performed in a restricted range of W and at higher Q^. A second fit restricted to the region 
1.5 < < 5 GeV^ at Q'^ = 2.4 GeV^, was also performed giving 6 = 0.44 ± 0.19(stat.). 
The fit is presented in Fig. [31 Also shown in the figure are previous ZEUS measurements at 
different values of [18] and the extension to higher W values from the present analysis. 
For each the corresponding S values fitted in the extended W range are given. Within 
the present accuracy the results do not show evidence for a dependence of 6. This 
result is similar to that obtained for the exclusive production of J/tp mesons [12-14,47]. 

The first direct measurement of the differential cross section da^''^'*"''^ / dt, extracted from 
the LPS-tagged events at = 3 2 GeV^ and at W = 104 GeV, is shown in Fig. H The 
value of the slope parameter b extracted from an exponential fit to the differential cross 
section, da^*P^^P/dt oc e"''!*!, is 6 = 4.5 ± 1.3(stat.) ± 0.4(syst.) GeV'^ (xVndf = 0.90). 
This value is consistent with the results obtained by HI [20] h = 5.45 ± 0.19(stat.) ± 
0.34(syst.) GeV~^ at = 8 GeV^ and = 82 GeV, from the transverse-momentum 
distribution of the photon candidate. 

A compilation of b values as measured for various exclusive processes [14,21], including 
the result of this paper, is shown in Fig. [5] as a function of + M^, where M is the mass 
of the exclusive final state. The b value presented here is lower but consistent with the 
corresponding vector mesons and HI DVCS values at similar scales. The fast rise of the 
DVCS cross section with W at Q"^ = 2.4 GeV^ and the low value of b at Q"^ = 3.2 GeV^ 
indicate that the DVCS process is a hard process even at low values. 

8 Summary 

The DVCS cross section has been measured as a function of and W in the region 
1.5 < < 100 GeV^ and 40 < IV < 170 GeV. The measured cross section decreases 
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steeply with Q^, showing a dependence Q"^", with n = 1.54 ± 0.05(stat.). The W cross 
section rises with increasing W following a functional form W^, with 6 = 0.52±0.09(stat.) 
and has little dependence on Q^. 

For the first time, the DVCS differential cross section as a function of t was measured by 
directly tagging the scattered proton. An exponential behaviour was assumed, the slope 
parameter b = 4.5 ± 1.3(stat.) ± 0.4(syst.) GeV~^ was obtained from a fit to the data at 
= 3.2 GeV^ and W = 104 GeV. These findings indicate that the DVCS process is a 
hard process even at low Q^. 
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g2 range (GeV^) 


g2 (GeV2) 


f^7*P-*7P (nb) 


1.5 - 5 


3.25 


21.28±0.92+}:°2 


5 - 10 


7.5 


5.87 ±0A2t',lt 


10 - 15 


12.5 


3.27 ±0.33^^6 


15 - 25 


20.0 


1.23 ±0.21^°:°^ 


25 - 40 


32.5 


0.55 ±0.18^|];J^^ 


40 - 100 


70.0 


0.16 ±0.07t|]:[]^ 



Table 1: The DVCS cross section, a"'''^^"'^, as a function of Q"^. Values are 
quoted at the centre of each bin and at W = 104 GeV . The first uncertainty is 
statistical and the second systematic. 





W range (GeV) 


W (GeV) 


^7*P-7P (nb) 


40 - 60 


50 


14.47±1.05+°J° 


60 - 80 


70 


20.38±1.57+}:°^ 


80 - 100 


90 


17.95±1.35l[J|^ 


100 - 120 


110 


20.65±1.26i;:?^ 


120 - 140 


130 


26.42±1.84+°]^ 


140 - 170 


155 


27.60±3.74+2;0i 



Table 2: The DVCS cross section, a'^ p^ip ^ ^ function of W . Values are 
quoted at the centre of each W bin and for = 3.2 GeV^ . The first uncertainty 
is statistical and the second systematic. 
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W range 


W 




^7*P-7P (nb) 


^7*P-7P (nb) 


^7*P-7P (nb) 


(GeV) 


(GeV) 


g2 = 2.4 GeV^ 


Q2 = 6.2 GeV^ 


Q2 = 9.9 GeV^ 


g2 = 18.0 GeV^ 


40 - 65 


52.5 


27.06±3.44tt42 








65 - 90 


77.5 


22.36±3.11+?:^° 








90 - 115 


102.5 


26.49±1.89;?:^^ 








115 - 140 


127.5 


35.94±2.63t}|^ 








140 - 170 


155 


35.72±9.47l|^^ 


16.93±2.43l{i^ 


6.15±1.67+[J:^} 


2.21±0.82l°:i^ 



Table 3: The DVCS cross section, cr^*P~''yp^ as a function ofW in four ranges. 
Values are quoted at the centre of each W bin and for the values listed. The 
first uncertainty is statistical and the second systematic. The values for higher 
and lower W , shown in Fig. are taken from a previous publication [18] and are 
not repeated here. 



da^*'P^^P/dt 


t range (GeV^) 


t (GeV2) 


a^'p-^ipjdt (nb/GeV^) 


0.08 - 0.19 


0.14 


34.6±9.6± 2.4 


0.19 - 0.31 


0.25 


32.7±9.4± 2.3 


0.31 - 0.42 


0.36 


19.6±7.5± 1.4 


0.42 - 0.53 


0.47 


5.7 ±4.1± 0.4 



Table 4: The DVCS differential cross section, da'^'^^'^P /dt, as a function of 
\t\. Values are quoted at the centre of each \t\ bin and for = 3.2 GeV^ and 
W = 104 GeV . The first uncertainty is statistical and the second systematic. 
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Figure 1: Distribution of (a) W , (h) in the inclusive sample and of (c) xl 
and (d) \t\ in the LPS sample, for the e-sample (top), the '^-sample (middle) and 
the '-y-sample after BH background and proton dissociation subtraction (bottom). 
Also shown are the expectations of the MC normalised to the luminosity of the data 
and the contribution from exclusive dilepton production (e'^e~ ). 
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Figure 2: (a) The DVCS cross section, u^'P'^^ip^ as a function of . The solid 
line is the result of a fit of the form ~ Q''^"' . (b) The DVCS cross section, a^'P^^P^ 
as a function of W . The solid line is the result of a fit of the form ~ . The 
inner error bars represent the statistical uncertainty while the outer error bars the 
statistical and systematic uncertainties added in quadrature. 
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Figure 3: The DVCS cross section, a"'*^'^'*^, as a function of W for = 
2.4 GeV^ (dots) shown together with previous ZEUS measurements (squares) [18]. 
Also shown at higher Q'^ are the new measurements at W = 155 GeV (dots). The 
solid lines are the results of a fit of the form a'^*P^'yp oc . The values of S and 
their statistical uncertainties are given in the figure. The inner error bars represent 
the statistical uncertainty while the outer error bars the statistical and systematic 
uncertainties added in quadrature. 
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Figure 4: The DVCS differential cross section, da"''^^"''^ /dt, as a function of 
\t\. The solid line is the result of a fit of the form ~ e~^'*'. The inner error bars 
represent the statistical uncertainty while the outer error bars the statistical and 
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Figure 5: A compilation of the values of the slope b as a function of Q"^ + 
for various exclusive processes including the present DVCS measurement. The 
inner error bars represent the statistical uncertainty while the outer error bars the 
statistical and systematic uncertainties added in quadrature. 
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